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Abstract 

The  student  under  the  training  program  is  first  taught  with  the  ground  work  of  crystal 
physics  and  the  transformation  of  axes  in  second-rank  tensors.  Then  he  has  participated  in 
different  topics  in  the  parent  research  grant  (F49620-92-J-0171),  "Mechanics  of  Localized 
Deformation  of  Solids  under  Cyclic  Mechanical  and  Thermal  Loadings". 

He  has  contributed  more  in  the  following  two  areas:  (1)  Effect  of  single  crystal  elastic 
anisotropy  on  the  fatigue  band  deformation  in  polycrystals.  The  application  of  Eshelby's 
Equivalent  Inclusion  Method  to  transform  an  inhomogeneous  solid  into  a  homogeneous  one  in 
shown  and  a  numerical  example  is  given.  (2)  The  analysis  of  the  growth  of  extrusions  in  a  single 
crystal  is  shown.  The  boundary  tractions  of  a  most  favorably  oriented  crystal  at  the  free  surface  of 
a  polycrystal  are  removed  to  simulate  a  single  crystal.  The  calculated  extrusions  are  compared  with 
those  observed  subgrain  boundary  movements. 
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I.  Introduction 


Single  crystals  under  loading  slide  along  certain  directions  on  certain  planes.  This  slip 
depends  on  the  resolved  shear  stress  in  the  slip  system  and  is  independent  of  the  normal  pressure 
on  the  sliding  plane.  This  crystal  characteristic  is  known  as  the  Schmid's  law.  This  Schmid's  law 
also  holds  under  cyclic  loading  (Parker,  1961).  In  a  face-centered  cubic  (fee)  material,  there  are 
four  slip  planes,  on  each  of  which  there  are  three  slip  directions,  giving  a  total  of  twelve  slip 
systems.  A  micromechanic  theory  of  fatigue  crack  initiation  has  been  developed  to  model  the 
gating  mechanism  of  fatigue  band  in  polycrystals  (Lin  and  Ito,  1969;  Lin  1992).  A  fatigue  band  is 
composed  of  three  thin  slices  P,  Q,  and  R  in  a  most  favorably  oriented  surface  crystal  (Figure  1) 
was  considered.  Lin  (1992)  has  shown  that  cyclic  tension  and  compression  loading  causes 
alternate  monotonic  buildup  of  plastic  shear  strains  in  both  P  and  Q  to  form  a  persistent  slip  band 
(PSB).  This,  in  turn,  causes  the  growth  of  extrusions  and  intrusions  (Forsyth  and  Stubbington, 
1955)  at  the  free  surface  of  the  polycrystal. 

These  PSB's  are  preferred  sites  of  fatigue  crack  initiation.  Hence,  the  depth  of  intrusion 
and  the  height  of  extrusion  are  taken  as  a  measure  of  the  fatigue  damage.  The  thickness  of  the 
slices  P,  Q,  euid  R  is  much  smaller  than  their  length  (along  X3-axis).  The  plastic  strain  caused  by 
slip  in  the  central  part  of  the  band  is  close  to  being  constant  along  the  length,  so  the  central  portion 
of  the  band  is  considered  to  be  under  plane  deformation. 

All  metals  have  initial  defects,  and  have  initial  stresses.  These  initial  stresses  play  an 
important  role  in  fatigue  crack  initiation.  Lin  and  Ito  (1969)  first  point  out;that  the  initial  stress 
field  favorable  to  extrusions  and  intrusions  can  be  caused  by  a  row  of  dislocation  dipoles.  Later, 
these  dislocation  dipoles  have  been  observed  in  the  ladder  structure  in  persistent  slip  bands 
(Antonopoulos,  1976;  Mughrabi,  1981).  These  dipoles  induce  initial  resolved  shear  stresses  at  the 
interface  between  the  PSB  and  the  matrix. 

Extrusion  causes  elongation  and  a  tensile  stress  in  R.  This  stress  combined  with  applied 
cyclic  stress  can  activate  a  second  slip  system  to  slide.  This  slip  in  the  secondary  slip  system 
enhances  the  extent  of  extrusion  (Lin  etal,  1989).  Interchanging  the  signs  of  the  initial  stresses  in 
P  and  Q  initiates  an  intrusion  instead  of  an  extrusion.  This  theory  of  fatigue  crack  initiation  has 
overwhelming  metallurgical  evidences  (Lin,  1977). 
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Grain 

Boundary 


Figure  1 

Fatigue  band  in  the  most  favorably  oriented  surface  crystal. 


11.  Elastic  Anisotropy  of  Crystal  in  Polycrystal 
Fatigue  Band  Formation 

The  above  model  mainly  concerned  with  the  fatigue  crack  initiation  of  aluminum  and  its 
alloys.  The  elastic  anisotropy  of  their  individual  crystals  is  insignificant  and  accordingly  neglected. 
Hence,  this  surface  crystal  is  assumed  to  be  embedded  in  an  elastically  isotropic  and  homogeneous 
semi-infinite  medium.  Thus  this  problem  is  reduced  to  an  anisotropic  inclusion  in  an  elastic 
isotropic  semi-infinite  solid.  This  heterogeneous  solid  is  then  transformed  into  a  homogeneous 
one  using  Eshelby's  Equivalent  Inclusion  Method  by  introducing  an  extra  set  of  eigenstrain 
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(Eshelby,  1957;  Mura,  1982).  Then  the  calculation  of  stress  and  strain  fields  can  be  carried  out  in 
the  same  way  as  for  a  single  phase  solid. 

A  method  of  inelastic  analysis  of  elastic  isotropic  homogeneous  solids  has  been  developed 
by  Lin  and  Ito  (1966)  and  Lin  (1971).  The  residual  stress  t^(x)  at  a  point  x  in  an  infinite  elastic 

isotropic  medium  caused  by  a  unit  plastic  strain  in  an  elemental  volume  Av'  at  x'  is 

expressed  as  (Lin,  1971) 


t:..(x)  =  GijTa{x,x')e^i(x')Av'  (1) 

where  Av'  is  an  elemental  volume  and  G^^(x,x')  is  the  influence  coefficient  of  the  stress  x~{x)  at 
X  caused  by  unit  plastic  strain  ^«(x')  in  x' .  In  the  following,  Av'  is  omitted  to  simplify  writing. 

Consider  first  a  homogeneous  isotropic  solid  subject  to  a  uniform  loading;  i.e.,  the  surface  traction 
Sf  along  Xj-axis  per  unit  area  with  an  exterior  normal  v  equals  Xg  Vj .  The  stress  is  uniform  and 

equal  to  Xg  with  a  corresponding  strain  eg .  If  this  body  has  subdomain  Q  in  the  region  D-Q,  as 
shown  in  Figure  2,  with  elastic  moduli  Cgj^i  different  from  those  of  D-Q  denoted  by  ^ijkl  >  this 
inclusion  would  give  a  disturbance  strain  eg  and  disturbance  stress  Xg .  Now  the  stress  in  Q 
would  be  Cgij{eg  +eg)  and  that  in  D-Q  would  be  Cg^ieg  +eg). 

Eshelby  (1957,  1961)  ingeniously  pointed  out  that  this  disturbance  stress  Xg  can  be 

* 

simulated  by  an  inelastic  strain,  referred  to  as  an  eigenstrain  eg  (Mura,  1982)  to  differentiate  it 
from  the  inelastic  strain  which  commonly  represent  the  plastic  and  creep  strain.  If  we  let 

'^ij  ~  Qjkl  (^kl  +  ^kl  ~^kl)~  ^ijkl  (^kl  ^kl )  ^ 

we  can  replace  the  subdomain  with  ^ijkl  by  one  with  Cg,j.  After  this  imaginary  replacement,  we 
have  a  homogeneous  solid  with  elastic  moduli  ^ijkl-  eg  has  the  same  effect  as  Cg  in  causing 
residual  stress.  We  have  from  Eq.  ( 1) 

Cyw[4(x)  -  ««(x)]  =  Gga(x,x')eli{x')  (3) 


Let  ^ijkl  be  the  elastic  compliance  of  the  region  D-Q.  From  Eqs.  (2)  and  (3), 


%(x)  =  Siju  Gga{^,x')eli(x’)  +  4(x') 


(4) 


Cijki4i  +  Cga[Skim„G„„rsi'^,x')ersix')  +  e«(x')J  =  +  Ggaix,x')euix')  (5) 
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Figure  2 

Anisotropic  surface  crystal  embedded  in  an  isotropic  semi-infinite  medium. 
This  gives  the  solution  of  %(x').  If  we  have  plastic  strain  e^,  Eqs.  (2)  and  (4)  become 


^ijkli^kl  ~^kl  +  ^kd  -  Qjkl  (^kl  “  +  ^kl ' 


4l) 


^ki('>0  =  Sijki  Gykiix,  x')[4(x')  +  4(x')]  +  4(x')  +  4(x') 
^jkl4l  +  +4(x')]+  4(x')} 

=  +  G^«(x>x')[4(x')  +4(x')] 


With  eg(x')  known,  Eq.  (8)  gives  the  solution  of  e^(x') . 


This  method  was  applied  to  find  the  plastic  strain  distribution  in  a  fatigue  band  in  the  most 
favorably  oriented  crystal  located  at  a  free  surface  of  the  polycrystal.  The  surface  crystal  was 

5(C 

considered  to  have  anisotropic  elastic  constants  and  the  remaining  part  of  the  polycrystal  was 
assumed  to  be  elastically  isotropic,  as  shown  in  Figure  2.  In  the  numerical  calculation,  the 
anisotropic  surface  crystal,  including  the  thin  slices  P,  Q,  and  R  was  divided  into  a  number  of 
elements.  Plastic  deformation  is  highly  concentrated  in  the  fatigue  band.  There  is  no  plastic  strain 
outside  P,  Q,  and  R.  By  using  the  method  described  above,  the  eigenstrain  e^j,  which  is  used  to 

transform  the  heterogeneous  body  into  a  homogeneous  one,  is  composed  of  two  parts,  one  from 
the  applied  load  and  the  other  from  the  plastic  strain.  The  contribution  of  these  parts  lead  to  the 
change  of  Schmidt  factors  and  stress  influence  coefficients,  and  hence  the  slip  distributions  in  the 
bands. 

In  the  following  numerical  calculation,  these  closely  located  thin  slices  P,  Q,  and  R  in  the 
crystal  are  considered  to  have  a  small  positive  initial  resolved  shear  stress  in  P  and  a  negative  one 
in  Q.  The  anisotropic  elastic  constants  of  NisAl  intermetallic  compound,  experimentally  obtained 
by  Yang  (1985),  are  used.  The  isotropic  elastic  constants  of  the  surrounding  metal  composed  of 
crystals  of  random  orientation  are  obtained  by  a  conventional  approximate  average  methods  such  as 
Voigt  or  Reuss  methods.  Assuming  the  strain-hardening  in  P  and  Q  to  be  zero  and  using  the 
numerical  method  developed  by  Lin  et  al.  ( 1989),  we  have  obtained  some  preliminary  plastic  strain 
distributions  in  the  fatigue  band  under  cyclic  tension  and  compression  at  10  and  20  cycles  as 
shown  in  Figure  3.  This  effect  of  elastic  anisotropy  will  be  further  studied  for  the  cases  of  multiple 
fatigue  band  in  the  most  favorably  oriented  crystal  and  on  the  propagation  of  the  fatigue  band 
across  the  grain  boundary.  The  heterogeneous  inclusion  causes  a  varying  stress  field  in  the  body 
under  uniform  loadings  X22  ■  The  Schmidt  factor  in  the  bands  is  0.5  if  the  solid  is  homogeneous 
and  isotropic.  With  the  anisotropic  inclusion  ^ijU  ,  the  Schimdt  factor  varies  along  the  bands.  The 

part  of  the  band  with  this  Schmidt  factor  greater  than  0.5  will  slide  and  start  to  form  a  fatigue  band. 
This  will  cause  different  distributions  of  plastic  strain  in  P  and  Q  and  affect  considerably  the 
growth  of  fatigue  band  of  the  material. 

The  microstress  of  the  fatigue  band  in  the  surface  crystal  depends  not  only  on  the 
anisotropic  elastic  constants  of  the  crystal,  but  also  on  these  anisotropic  constants  of  all  the 
crystals.  The  previous  study  considers  only  the  elastic  anisotropy  of  the  surface  crystal.  Further 
study  considers  also  the  elastic  anisotropy  of  the  surrounding  crystals.  The  details  of  the  analysis 
and  the  results  are  given  in  the  paper  "Elastic  Anisotropy  Effect  of  Crystals  on  Polycrystal  Fatigue 
Crack  Initiation",  published  in  the  Journal  of  Engineering  Materials  and  Technology,  ASME, 
Vol.  1 17,  No.  4,  pp.  470-477.  A  copy  of  the  paper  is  given  in  the  appendix. 
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Figure  3 

Plastic  shear  strain  distribution  along  band  considering  band  anisotropy  in  the  surface  crystal. 


III.  Fatigue  Bands  in  Single  Crystal 

Single  crystals  are  often  used  in  jet  engines.  Extrusions  and  intrusions  in  PSBs  have  been 
observed  in  these  single  crystals  under  high-cycle  fatigue  (HCF)  loadings.  These  PSBs  are  the 
preferred  sites  of  fatigue  crack  initiation.  Hence  this  study  is  of  both  scientific  interest  and 
engineering  needs. 

The  growth  of  extrusion  and  intrusion  in  polycrystal  of  elastic  isotropy  has  been  studied  by 
Lin  and  his  associates  (Lin  etal.,  1989;  Lin,  1992).  To  analyze  this  growth  in  a  polycrystal,  the 
fatigue  bands  are  assumed  to  occur  in  a  most  favorably  oriented  crystal  embedded  at  a  free  surface 
of  the  metal.  The  analysis  is  reduced  to  a  plane  strain  solution  of  a  semi-infinite  elastically 
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isotropic  solid.  The  cross-section  of  the  single  crystal  is  assumed  to  be  rectangular  (Figure  4)  with 
"b"  considerably  larger  than  "a".  The  deformation  of  the  center  portion  will  be  close  to  plane 
strain.  The  fatigue  band  in  this  portion  is  considered  to  represent  that  of  the  single  crystal.  Using 
the  solution  of  fatigue  bands  in  polycrystal  (Lin,  1992),  we  obtain  the  stress  field  of  this  surface 
crystal  with  boundary  tractions  as  shown  in  Figure  5(a).  To  find  the  stress  field  due  to  this  slid 
grid  in  the  single  crystal,  these  boundary  tractions  are  to  be  relieved  by  a  set  of  equal  and  opposite 
tractions  using  finite  element  method  (FEM),  as  shown  in  Figure  5(b).  The  sum  of  the  stresses 
given  by  Figures  5(a)  and  5(b)  gives  the  residual  stress  field  due  to  the  slid  grid,  as  shown  in 
Figure  5(c). 
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Figure  4 

Single  crystal  under  cyclic  tension  and  compression  loadings. 


(a)  (b)  (c) 

Figure  5 

Superposition  of  semi-infinite  plane  strain  solution  and  finite  element  solution. 
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The  localization  of  the  surface  traction  (Figure  5(b))  depends  on  the  distance  of  the  slid  grid 
from  the  boundary.  The  amount  of  numerical  work  and  the  difficulty  of  the  FEM  to  achieve  a 
given  accuracy  depends  on  the  distance.  Hence,  for  a  slid  grid  in  the  left  half  of  the  crystal,  as 
shown  in  Figure  5(a),  the  free  surface  of  the  semi-infinite  solid  is  taken  on  the  left.  Similarly,  for  a 
slid  grid  in  the  right  half,  the  free  surface  is  taken  to  be  on  the  right.  By  this  way,  the  boundary 
traction  would  be  much  less  localized.  This  would  facilitate  greatly  the  FEM  calculations. 

In  1980,  Mecke  and  Blochwitz  performed  an  experiment  using  nickel  single  crystal  to 
study  the  behavior  of  PSBs  under  cyclic  loadings.  This  experiment  was  conducted  at  room 
temperature  and  at  a  constant  plastic  strain  amplitude  of  =  l.SxlO'^.  The  result  is  shown  in 

Figure  6  after  many  cycles  of  loadings.  This  figure  shows  that  both  the  PSBs  and  the  subgrain 
boundaries  (SGB),  1  to  2  p,m  each  in  width,  have  penetrated  across  the  whole  crystal  and  extruded 
from  both  sides. 


P5B  PSB  PSB 


Using  the  numerical  method  developed  as  described,  we  consider  a  seven  band  face- 
centered  cubic  single  crystal  loaded  in  cyclic  tension  and  compression.  This  model  is  shown  in 
Figure  7.  The  most  favorably  slip  system,  called  the  primary  slip  system,  makes  an  angle  of  45° 
with  the  specimen  axis.  As  in  the  previous  models  (Figure  1),  each  fatigue  band  is  composed  of 
three  thin  slices.  These  slices  in  the  left  half  are  denoted  by  P,  Q,  and  R,  and  those  in  the  right  half 
are  denoted  by  P',  Q',  and  R'.  As  R  and  R'  extrudes,  a  second  slip  system  may  be  activated. 
This  crystal  is  assumed  to  be  elastically  isotropic.  The  shear  modulus  p  is  taken  to  be  50  GPa  and 
Poisson's  ratio  is  taken  to  be  0.3.  The  critical  shear  strength  0“^  is  assumed  to  be  200  MPa  and  the 
applied  stress  022  is  400  MPa. 
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Figure  7 

Isotropic  seven  band  model  of  the  single  crystal. 

To  see  the  effect  of  the  location  of  initial  tensile  strain  in  R  on  the  extrusion  growth,  the 
grid  with  this  initial  strain  was  taken  to  be  in  the  left,  then  at  the  center,  and  finally  at  the  right 
portion  of  a  single  fatigue  band.  The  calculated  plastic  strain  variations  with  the  number  of  loading 
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cycles  are  shown  in  Figure  8.  All  three  locations  show  positive  shear  strain  in  P  and  a  negative 
one  in  Q,  showing  the  occurrence  of  extrusions. 


Figure  8 

Plastic  strain  distribution  at  left  free  surface. 

Then  the  initial  tensile  stress  at  the  center  and  the  applied  cyclic  loadings  are  increased  until 
the  plastic  strain  occur  in  both  the  primary  and  secondary  slip  systems.  These  primary  and 
secondary  plastic  strains  variations  along  the  band  at  different  loading  cycles  are  shown  in  Figure  9 
and  10,  respectively. 
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Figure  9 

Plastic  strain  in  primary  slip  system  of  Q. 
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Figure  10 

Plastic  strain  in  secondary  slip  system  of  R. 


These  gives  some  preliminary  results  of  HCF  bands  in  single  crystal. 
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Fatigue  bands  have  been  observed  in  both  monocrystalline  and  polycrystalline  metals. 
Extrusions  and  intrusions  at  the  free  surface  of  fatigued  specimens  are  favorable 
sites  for  fatigue  crack  nucleation.  Previous  studies  (Lin  and  Ito,  1969;  Lin^  1992) 
mainly  concerned  the  fatigue  crack  initiation  in  aluminum  and  its  alloys.  The  elastic 
anisotropy  of  individual  crystals  of  these  metals  is  insignificant  and  was  accordingly 
neglected.  However,  the  anisotropy  of  the  elastic  constants  of  some  other  metallic 
crystals,  such  as  titanium  and  some  intermetallic  compounds,  is  not  negligible.  In 
this  paper,  the  effect  of  crystal  anisotropy  is  considered  by  using  Eshelby's  equivalent 
inclusion  method.  The  polycrystal  analyzed  is  NUAl  intermetallic  compound.  The 
plastic  shear  strain  distributions  and  the  cumulative  surface  plastic  strain  in  the 
fatigue  band  versus  the  number  of  loading  cycles  were  calculated,  and  the  effect  of 
crystal  anisotropy  on  the  growth  of  the  extrusions  was  examined. 


1  Introduction 

The  process  of  fatigue  failure  in  ductile  metals  is  divided 
into  the  following  three  stages:  ( 1 )  the  crack  initiation,  (2)  the 
growth  and  coalescence  of  rhicrocracks  to  form  some  mac¬ 
rocracks,  and  (3)  the  propagation  of  dominant  cracks  leading 
to  the  final  catastrophic  failure.  In  high-cycle  fatigue,  fatigue 
crack  initiation  predominates  fatigue  life.  Fatigue  bands  are  the 
favorable  sites  for  the  nucleation  of  fatigue  microcracks. 

Fatigue  bands  consist  of  highly  localized  plastic  strains  and 
have  been  observed  in  both  single  crystals  and  polycrystals. 
Plastic  strains  are  mainly  due  to  the  movement  of  dislocations. 
Under  the  cyclic  loading,  the  localized  plastic  deformation  in¬ 
creases  with  the  number  of  loading  cycles  and  causes  extrusions 
and  intrusions  at  the  free  surface  of  metals.  The  height  of  extru¬ 
sions  or  the  depth  of  intrusions  is  here  taken  as  a  measure  of 
fatigue  damage. 

A  theory  of  micromechanics  of  fatigue  crack  initiation  has 
been  proposed  to  model  the  gating  mechanism  of  fatigue  bands 
in  polycrystals  (Lin  and  Ito,  1969;  Lin,  1992).  A  fatigue  band 
composed  of  three  thin  slices  P,Q,  and  Rim  most  favorably 
oriented  surface  crystal  was  considered.  Under  cyclic  tension 
and  compression  loading,  the  alternate  monotonic  build-up  of 
plastic  shear  strains  in  both  P  and  Q  causes  the  growth  of 
extrusions  or  intrusions  at  the  free  surface  of  polycrystals.  Since 
this  model  is  concerned  mainly  with  the  fatigue  crack  initiation 
of  aluminum  and  its  alloys,  the  elastic  anisotropy  of  their  indi¬ 
vidual  crystals  is  insignificant  and  accordingly  neglected.  Hence 
the  surface  crystal  is  assumed  to  be  embedded  in  an  elastically 
isotropic  and  homogeneous  semi-infinite  medium.  For  some 
other  polycrystalline  metals,  such  as  titanium  and  some  inter- 
metalic  superalloys,  however,  the  elastic  anisotropy  of  crystals 
is  not  negligible.  For  such  polycrystals,  the  local  microstress 
field  in  a  crystal  depends  not  only  on  the  applied  stresses  but 
also  on  the  orientations  of  its  surrounding  crystals.  Hence  the 
elastic  anisotropy  of  crystals  will  affect  the  fatigue  crack  initia¬ 
tion  in  these  metals. 

In  this  paper,  the  elastic  anisotropy  of  individual  grains  in  a 
polycrystal  is  considered.  A  fatigue  band  is  assumed  to  occur 
in  ihe  most  favorably  oriented  surface  crystal.  This  crystal  is 
surrounded  by  the  crystals  of  various  given  orientations.  The 
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crystals  outside  the  surrounding  crystals  are  randomly  oriented 
and  assumed  to  be  isotropic  and  homogeneous.  Hence  this  sur¬ 
face  crystal  and  the  surrounding  crystals  of  given  orientations 
are  considered  to  be  embedded  in  a  semi-infinite  isotropic  and 
homogeneous  elastic  medium.  Thus  it  becomes  an  inclusion 
problem.  This  heterogeneous  solid  is  transformed  into  a  homo¬ 
geneous  one  by  using  Eshelby’s  equivalent  inclusion  method 
(Eshelby,  1957,  1961;  Mura,  1982).  In  this  method,  an  extra 
set  of  eigenstrains  is  introduced.  These  eigenstrains  can  be  ex¬ 
plicitly  expressed  in  terms  of  the  applied  stresses  and  plastic 
strains.  After  this  transformation,  the  procedures  for  this  micro¬ 
mechanical  analysis  are  the  same  as  those  given  by  Lin  ( 1992) 
for  elastically  isotropic  polycrystals.  In  the  following  numerical 
calculation,  ihe  anisotropic  elastic  constants  of  NisAl  monociys- 
talline  intermetallic  compound  (Yang,  1985)  were  used.  The 
isotropic  elastic  constants  of  the  surrounding  metal  were  taken 
from  the  experimental  data  (Stoloff,  1989).  The  plastic  strain 
distributions  and  the  cumulative  surface  plastic  strains  in  the 
fatigue  band  versus  the  number  of  loading  cycles  were  calcu¬ 
lated,  and  the  effect  of  elastic  anisotropy  on  the  formation  of 
extrusions  and  intrusions  was  examined. 

2  Gating  Mechanism  Provided  by  Microstress  Field 

Single  crystal  tests  show  that  slip  occurs  on  certain  crystallo¬ 
graphic  planes  along  certain  directions.  This  slip  depends  on 
the  resolved  shear  stress  in  the  slip  system  and  is  independent 
of  the  normal  stress  on  the  sliding  plane.  The  dependence  of 
slip  on  the  resolved  shear  stress  is  known  as  the  Schmid  law. 
It  has  been  shown  that  the  Schmid  law  also  holds  when  crystal¬ 
line  solids  are  subjected  to  cyclic  loading  (Parker,  1961).  For 
a  face-centered  cubic  ciystal  (FCC),  there  are  four  slip  planes, 
on  each  of  which,  there  are  three  slip  systems.  These  twelve 
slip  systems  are  shown  in  Fig.  1. 

All  metals  have  initial  defects,  hence,  have  initial  stresses. 
These  initial  stresses  play  an  important  role  in  fatigue  crack 
initiation.  Lin  and  Ito  ( 1969)  first  point  out  that  the  initial  stress 
field  favorable  to  extrusions  and  intrusions  can  be  caused  by  a 
row  of  dislocation  dipoles.  Later,  these  dislocation  dipoles  have 
been  observed  in  the  ladder  structure  in  a  persistent  slip  band 
(PSB)  (Antonopoulos,  1976;  Mughrabi,  1982).  These  dipoles 
induce  initial  resolved  shear  stresses  at  the  interface  between 
the  PSB  and  the  matrix. 

Based  on  experimental  observations,  a  micromechanic  theory 
of  fatigue  crack  initiation  was  proposed  by  Lin  and  Ito  (1969). 
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Fig.  1  Crystallographic  slip  planes  and  directions  of  FCC  crystals 


In  this  model,  a  most-favorably  oriented  surface  crystal  con¬ 
taining  a  fatigue  band  composed  of  three  thin  slices  P,  Q,  and 
/?  in  a  semi-infinite  poly  crystal  shown  in  Fig.  2  is  considered. 
The  slip  plane  and  slip  direction  of  the  primary  slip  system  a/? 
form  an  angle  of  45  deg  with  the  loading  axis.  Lin  and  Lin 
(1983)  have  calculated  a  uniform  initial  resolved  shear  stress 
field  positive  in  P  and  negative  in  Q,  caused  by  a  distribution 
of  plastic  strains.  This  initial  stress  field  satisfies  the  conditions 
of  both  compatibility  and  equilibrium.  During  the  cyclic  load¬ 
ing,  plastic  strains  occur  in  P,  Q  and  R  and  cause  a  residual 
stress  field.  The  total  resolved  shear  stress  r  in  a  slip  system  is 
the  sum  of  three  stresses  as 


7-  =  -1- 

where  is  the  resolved  shear  stress  caused  by  the  applied  load 
and  the  residual  resolved  shear  stress  due  to  plastic  strains. 

A  tensile  loading  ah  on  the  polycrystal  produces  a  positive 
in  the  entire  surface  crystal.  Since  rMn  P  is  positive,  the 
total  resolved  shear  stress  r  in  P,  which  is  reaches 

the  critical  shear  stress  first;  and  hence,  according  to  Schmid 
law,  P  slides.  The  plastic  shear  strain  caused  by  this  slip  in  P 
produces  a  residual  stress  field  which  has  negative  resolved 
shear  stress  in  Q  and  makes  Q  to  slide  more  readily  in  the 
reverse  loading.  The  negative  slip  in  Q  induces  in  turn  a  positive 
residual  resolved  shear  stress  in  P  making  P  more  readily 
to  slide  under  the  next  tensile  loading.  This  process  is  repeated 
during  the  cyclic  loading  for  every  cycle  thus  providing  a  natural 
gating  mechanism  for  a  monotonic  buildup  of  local  resolved 
plastic  shear  strain  in  both  P  and  Q,  pushing  R  out  of  the 
free  surface  and  promoting  an  extrusion.  The  elongation  of  R 
may  activate  a  secondary  slip  system  to  slide.  The  effect  of  slip 
in  the  secondary  slip  system  enhances  the  extent  of  extrusion 
(Lin,  et  al.,  1989)  .  Interchanging  the  signs  of  the  initial  stresses 
in  P  and  Q  initiates  an  intrusion  instead  of  an  extrusion.  This 
micromechanical  theory  on  fatigue  crack  initiation  has  over¬ 
whelming  metallurgical  experimental  evidences  (Lin,  1977). 

3  Calculation  of  Residual  Stress  Field 

3.1  Residual  Stress  Field.  The  residual  stress  field  due 
to  a  distribution  of  inelastic  strains  can  be  calculated  by  the 
analogy  between  inelastic  strains  and  applied  forces  (Lin, 
1969).  It  has  been  shown  that  the  equivalent  body  force  per 
unit  volume  along  the  x^-axis  due  to  inelastic  strains  elj  is 

Fi  =  -C^kiekj  (1) 

where  Q*/  is  the  elastic  modulus.  The  summation  convention 
is  used  here  and  the  subscript  after  a  comma  denotes  the  differ¬ 
entiation  with  respect  to  the  coordinate  variable.  For  isotropic 
elastic  solids,  the  elastic  constants  are 

Cijki  =  ^SijSki  +  fj,(6ikSji  +  SiiSjk) 

where  \  and  //  are  Lame’s  constants  and  Sij  is  the  Kronecker 
delta.  The  equivalent  surface  force  per  unit  area  along  the  x,- 
axis  has  been  shown  as 


n  j  I  n  1 1 1 


ii  ii  I  n  n 


Si  —  Cijkie^Vj 


(2) 


where  vj  is  the  direction  cosine  of  the  normal  to  the  area.  The 
fictitious  stresses  o’J(x)  at  point  x  caused  by  equivalent  forces 
in  domain  H  can  be  written  by  using  the  Green’s  function  as 

=  f  T'tj(.x,x')Ft(x')dQ,  +  f  T^(x,  x')St(x')dr 

Jn  Jr 

=  -Ca™.  rj(x,  x')€'^j(x')dQ 

+  r^j(x,  x')e'^(x')i'!dr  (3) 


where  r|(x,  x')  represents  the  stress  component  at  point  x 
induced  by  a  unit  force  acting  at  point  x'  along  the  x*-axis.  T 
is  the  boundary  of  domain  Q.  If  inelastic  strain  components  of 
e^(x)  are  constant  in  O,  eqn  (3)  is  reduced  to 


aijix)  =  tJ(x,  x')Mrj€:„(j^)  (4) 


where  e  "„(n)  is  the  constant  inelastic  strain  in  Cl.  The  residual 
stress  field  is 


Fig.  2  Fatigue  band  in  the  most  favorably  oriented  surface  crystal 


ctJ(x)  =  a1j(x)  -  Q*,eS(x)  (5) 
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Fig.  3  Euier  angles  and  representation  of  crystal  orientations 


3.2  Residual  Stress  Influence  Coefficient.  To  calculate 
the  residual  stress  field  numerically,  the  discretization  of  a  con¬ 
tinuous  problem  is  needed.  In  the  following,  the  domain  O  in 
an  extended  solid  is  divided  into  a  number  of  subdomains  or 
elements  In  each  of  them,  it  is  assumed  that  the  inelastic 
strain  e,"  is  uniform.  After  applying  Eqs.  (4)  and  (5)  to  the 
residual  stress  field  due  to  €jJ,:(x')(x'  €  0^)  is 


CTijW 


|J^  C^TUx,x')u,dr,  -  (6) 


where  the  function  JT(x,  is  defined  as 


H(x,  n,) 


1,  X  G  0^ 
0,  X  ^  n. 


Generally,  residual  stress  field  a-J(x)  varies  from  one  point  to 
another.  The  volume  average  of  aj(x)  over  subdomain  Clp  is 
taken  to  represent  the  residual  stresses  in  element  Qp.  Letting 
cT^jiCtp)  represent  the  residual  stresses  in  flp  due  to  in 
we  have 


<^22 

}  }  »  )  I 


>  >  <  >  } 

<^22 


cr^jiQp)  =  Gymninp.  (7) 

where  Gijmni^p,  ^9)  are  called  the  residual  stress  influence  coef¬ 
ficients,  which  represent  the  residual  stresses  in  subdomain  Qp 
caused  by  unit  inelastic  strain  components  in  subdomain 
For  the  resolved  residual  shear  stress  in  slip  system  afi 
in  Clp  due  to  the  resolved  plastic  shear  strain  e^rj  in  slip  system 
^77  in  n,,  Eq.  (7)  can  be  rewritten  by  the  tensor  transformation 
as 


K^i^p)  =  Gia/3,  Qp;  ^77,  (8) 

where  G(ap,  represents  the  residual  resolved  shear 

stress  in  slip  system  ap  in  caused  by  a  unit  resolved  plastic 
shear  strain  in  slip  system  ^77  in  The  total  resolved  shear 
stress  Tafi  in  slip  system  aP  in  is 

TaB{%)  =  ^'^0  +  +  X  X  Q.p\ 

q  in 

(9) 

The  summation  with  regard  to  slip  system  ^77  includes  all  the 
slip  systems  in 

3.3  Generalized  Plane  Strain  Problem.  The  thickness 
of  fatigue  bands  is  much  smaller  than  the  length  (dimension 
along  the  jcs-axis)  as  observed  at  the  free  surface  of  metals. 
Therefore,  the  problem  is  considered  under  the  plane  strain 
deformation.  The  plastic  strain  due  to  slip  in  the  secondary  slip 
system  is  not  confined  to  the  X1X2  plane.  This  causes  an  equiva¬ 
lent  force  corqponent  F3  acting  along  the  xs-direction.  The  exis¬ 
tence  of  this  F3  requires  a  modification  of  the  plane  strain  prob¬ 
lem.  Hence,  the  generalized  plane  strain  problem  is  considered 
in  the  present  study.  A  similar  problem  is  shown  in  the  analysis 
of  prismatic  anisotropic  bars  by  Lekhnitski  ( 1963).  This  gener¬ 
alized  plane  strain  deformation  is  defined  as 

Ui  —  Ui(xu  X2),  i  =  1,  2,  3 

where  Ui  is  the  displacement  component  along  X/-axis.  This 
gives 

Oij  =  2^1^  6ij6  +  i  (Uij  +  Ujj)^ 


Fig.  4  Polycrystal  model  of  fatigue  crack  initiation 


where  ^  =  mi,i  +  M2.2  and  u  is  the  Poisson’s  ratio.  The  equilib- 
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Fig.  6  Modified  Schmid  factor  varying  along  slice  P  Fig.  8  Plastic  shear  strain  distributions  in  P  at  100  cycles  and  2000 

cycles 


rium  equations  in  terms  of  displacement  components  m,  can  be 
expressed  as 


and 


v\- 


1  de  ^ 

\  —  2(j.  dxi  jx 


0,  i  =  1,  2 


(10) 


and 


ru(x,  x')  =  t^3(x,  x')  =  0,  *=1,2  (12) 


rfsCx,  x')  = 


dxi 


r|3(x,x')  = 


dx2 


where 


V2„3  +  =  0 


(11) 


dx\  dX2 


and  Fi  is  the  body  force  component.  Equations  (10)  and  (11) 
are  not  coupled  and  can  be  solved  separately. 

The  stress  components  r”(x,  x')  in  Eq,  (3)  for  the  semi¬ 
infinite  medium  with  the  generalized  plane  strain  deformation 
have  been  expressed  in  terms  of  the  Airy’s  stress  functions  (Lin 
and  Lin,  1974)  as 


dxl 


r?,(x,x')  =  l^,  rkx,x')  =  "'"^ 


dx\  ’ 


rt2(x,  X')  =  -  ,t^3(x,  x')  =  vS/^<i>k, 

UX\  UX2 


r|(x,  x')  ^0,  for  other  i  and  j  (13) 

where  the  AJry’s  stress  functions  are 
</)i(x,  x')  =  “(p  +  q)(x2  -  X2)(0i  +  02) 

+\q{x,  -  x\)  In  (X./Xz)  +  2px,x\{xx  +  x[)IX2 
<^2(x,  x')  =  ip  +  q)(x2  -  x|)(i9,  +  ^2) 

+  \q{x2  -  x'2)  In  (X./Xj)  -  2px,x\{x2  -  ;cD/X2 
(^3(x,  x')  =  -  —  (In  X,  +  In  X2) 


P  = 


1 


47r(l  -  i/) 


61  =  arctan 


^2  -  ^2 


q^p{l  -  2v) 


-TT  <  <  TT 


Fig.  7  Plastic  shear  strain  at  free  surface  versus  number  of  cycles 


Xi  =  (Xi  -  xiy  +  (X2  -  xi)^, 
X2  =  (Xi  +  x[)^  +  (,X2  -  xi)^ 


By  using  Eqs.  (6),  (12)  and  (13),  the  residual  stresses  as 
well  as  the  residual  stress  influence  coefficients  due  to  the  in¬ 
elastic  strains  in  the  element  of  a  specific  shape  can  be  readily 
calculated. 

4  Generalization  of  Equivalent  Inclusion  Method 

Eshelby’s  equivalent  inclusion  method  (Eshelby,  1957, 
1961 )  is  here  used  to  transform  an  inhomogeneous  problem  into 
a  homogeneous  one  by  introducing  a  distribution  of  eigenstrains 
(Mura,  1982).  Originally,  this  method  was  applied  to  solve  the 
problems  with  an  ellipsoidal  inhomogeneity.  In  the  present 
study,  this  method  is  generalized  to  solve  inhomogeneous  prob¬ 
lems  of  arbitrary  shapes. 
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In  the  following,  the  matrix  and  vector  notation  will  be  used 
instead  of  tensor  notation  for  convenience.  It  is  easy  to  contract 
a  second-order  tensor  with  components  Gy*/  into  a  matrix  with 
entries  G/,-  by  setting,  for  instance,  Gn  =  Gnu,  Gu  =  Gnizy 
Oi6  =  <^1112  +  Gu2i  and  G^i  =  Gmu  etc.  Vectors  are  denoted 
by  boldface  lowercase  Roman  or  Greek  letters  and  matrices  by 
boldface  uppercase  Roman  letters.  The  Eshelby’s  equivalent 
equation  is 

C*(x)[€°  +  e(x)  ~  e^(x)] 

=  C[€°  +  €(x)  -  €*(x)  -  e^(x)],  X  €  n  (14) 

where  H  is  an  inclusion  or  a  subdomain  with  the  elastic  constant 
C*(x)  embedded  in  an  extended  body  D  with  the  isotropic 
elastic  constant  C.  €(x)  is  the  disturbance  strain,  €*(x)  the 
eigenstrain,  e^(x)  the  plastic  strain  and  e®  uniform  strain  from 
which  =  Cc®  where  o'®  is  the  applied  far-field  stress.  The 
stress  O'  and  strain  e  are  vectors  expressed  as 

^  0‘22j  O'ss,  0'23,  0'i3,  0'i2}^ 

^  ~  {^II,  e22>  ^33,  €23»  €12}^- 

where  notation  {  •  •  • }  ^  represents  the  transposition.  Equation 
(14)  shows  that  the  eigenstrain  €*(x)  plays  the  same  role  as 
the  plastic  strain  e^(x)  does  in  the  homogeneous  materials. 

When  a  continuous  problem  is  discretized  for  the  purpose  of 
numerical  analysis,  the  domain  of  Q  is  further  divided  into  N 
subdomains  n,  ( 1  <  f  <  N)  called  elements,  in  each  of  which 
both  the  eigenstrain  e*  and  plastic  strain  are  assumed  to  be 
uniform.  The  stresses  in  element  Cli  caused  by  the  unit  eigens¬ 
train  €*  in  element  Qj  are  the  same  as  those  caused  by  the  unit 
plastic  strain  Thus  the  stress  in  induced  by  both  of 
the  eigenstrains  and  plastic  strains  in  all  the  elements  is 

o-f  =  X  G,ye;  +  I  l^i^N  (15) 

J  k 

where 


*_r*  *  *  *  *  *ir 

—  i^llj  ^22^  €33,  e23,  ^12); 


P  f  P  P  P  P  P 
€j  —  {Cii,  €22*  €33,  e23>  ^13,  € 


P  iT 
12  }j 


in  Qj  and  Gy  are  6  X  6  matrices  of  the  residual  stress  influence 
coefficients,  which  are  contracted  from  G*/^(n,  ,  Q^). 

For  the  generalized  plane  strain  problem,  W3  does  not  vary 
along  the  Xs-direction.  The  total  normal  strain  component  is 
equal  to  zero.  The  strain  component  €33  is  here  taken  to  be  zero. 


6.0- 

5.0- 

4.0- 

3.0- 


z^^±DAmPa 
=  IQOMPa 


Isotropic  Cose 


Anisotropic  Cose 


200  300 

Number  of  Cycles 


Fig,  9  Comparison  of  surface  plastic  strains  with  and  without  consider¬ 
ing  crystal  anisotropy 


Since  the  total  strain  €33  is  the  sum  of  £33  and  £33,  hence  £33  = 
0.  The  stress-strain  relation  may  be  given  as 

Ce,  -  C(€f  -h£f)  (16) 

where 

X  +  2fx  X  0  0  0 

X  X  -I-  2/i  0  0  0 

X  X  0  0  0 

0  0  2//  0  0  ’ 

0  0  0  2jji  0 

-  0  0  0  0  2^. 

X  +  2/i  X  X  000“ 

X  X  +  2/i  X  000 
X  X  X  +  2/i  0  0  0 

0  0  0  2/i  0  0 

0  0  0  0  2/i  0 

0  0  0  0  0  2/i . 


c  - 

and 


=  {en,  £22,  ^23,  ei3,  Cu)!  with  £33  ^  0 


From  Eqs.  (15)  and  (16),  expressing  the  disturbance  strains  €, 
in  Qi  in  terms  of  ef  and  €*,  we  obtain 


=  C-^[C(6*  +  €?)  +  I  Gye/  +  I  G,*€f], 


where 


1  <  i  <  V  (17) 


1  -1/  -u  0  0  0 

-1/  I  -1/  0  0  0 

0  0  0  1  +  0  0 

0  0  0  0  1  +  z/  0 

0  0  0  0  0  1  +  i/_ 

in  which  E  is  the  Young’s  modulus. 

The  left-hand  side  of  Eq.  (14)  can  also  be  written  for  the 
generalized  plane  strain  problems  as 


cr,  =  o-f  +  C*€,-  -  C*ef  (18) 

where  C*  is  a  6  X  6  anisotropic  elastic  stiffness  matrix  and  C* 
a  6  X  5  matrix  formed  from  matrix  C*  by  taking  the  third 
column  away.  And  erf  =  C*e®  in  which 

{^ii,  £22»  ^23,  Cl3»  €12}^  with  €33  ^  0 

Substituting  Eq.  (17)  into  Eq.  (18)  and  Eq.  (15),  Eq.  (14) 
becomes 


T€*  +  I  S,€;  =  q. 


+  (C*-T)€f-IS,*€j, 

k 


lsi:£N  (19) 


where 


T  =  C*C"'C 
Sy  =  (C*C-'-I)Gy 


and  I  is  a  6  X  6  unit  matrix. 

It  can  be  seen  that  eigenstrains  ef  depend  on  the  applied 
load  O'®  and  the  distribution  of  plastic  strains  €*  as  well.  Equa¬ 
tion  (19)  is  a  set  of  linear  algebraic  equations  with  unknowns 
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ef .  Solving  these  equations  for  ef  in  all  elements  ,  we  can 
express  ef  as 

ef  =  p.-  +  I  l^i^N  (20) 

* 

where  p,  is  a  vector  in  terms  of  cr°  in  and  Rik  is  the  matrix 
relating  the  plastic  strains  ei  in  to  the  eigenstrains  e  f  in  . 
The  total  stress  in  Q,  ,  obtained  by  substituting  Eq.  (20)  into 
Eq.  (15),  is 

o-i  =  +  O’- 

=  cr®  +  X  ^yP;  +  X  (^iit  +  S  GijR.jk)€{  (21) 
j  *  J 

In  the  above  equation,  the  term  of  2  GyP;  is  a  modification 
to  the  effect  of  the  applied  load  and  2  GyR/*  the  one  to  the 
matrix  G,*  of  the  residual  stress  influence  coefficients  due  to 
the  inhomogeneity. 

Under  the  stress  component  0-22  alone,  Eq.  (21)  reduces  to 
o-i  =  P,  cr22  +  X  G,jt€* 

k 

where 

p.(J^  =  £7  +  2](  GyP; 

J 

and 

Giic  =  Gik  +  X  GyR;fc 

J 

For  slip  system  a/3  in  element  fli ,  we  have 

€i  —  Olfxp€a0i 


and 

~  ^apOj 

where  ej/?,  and  are  the  plastic  shear  strain  and  resolved 
shear  stress,  respectively,  and 

T^aP  =  ^22^  "^33,  ^23,  ^13^  'Wll  1 1/3 

^ap  ~  {Wll>  ^22^  ^33 ^  ^23^  ^13*  ^l2}aP 

in  which  wiy  and  riij  are  composed  of  direction  cosines  of  slip 
system  a/B.  The  total  resolved  shear  stress  in  the  slip  system  is 

T f^Pi  “  aPi  api  api  (22) 

in  which 

'T'a/Jj  =  T^^pPi^22  ~  Sapi<^22 

^  X  X  G(a/?,  Q>j\  ^7], 

k  ey 

where  s^pi  is  the  modified  Schmid  factor  with  the  consideration 
of  anisotropic  elastic  constants  and  G(aP,  f2/ ;  ^77,  f^;)  expressed 
as 

G(apj  ;  ^77,  Q*)  =  n^pGikm^rj 

are  the  modified  residual  stress  influence  coefficients  after  trans- 
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forming  the  inhomogeneous  problem  into  a  homogeneous  one. 
The  summation  with  regard  to  slip  system  ???  includes  all  the 
slip  systems  in  /:th  element.  It  can  be  seen  that  Eq.  (22)  corre¬ 
sponds  to  Eq.  (9). 

5  Effect  of  Crystal  Elastic  Anisotropy  on  Fatigue 
Crack  Initiation  in  Polycrystals 

5.1  Elastic  Constants  of  Monocrystalline  and  Polycrys¬ 
talline  Crystals.  In  the  present  study,  nickel-base  superalloy 
polycrystals  of  FCC  structure  are  taken  for  the  numerical  analy¬ 
sis.  Fatigue  crack  initiation  of  other  materials  can  be  analyzed 
by  the  same  approach. 

There  are  many  experimental  data  available  for  nickel-base 
superalloy  single  crystals  in  the  literature.  Since  the  alloy  com¬ 
positions  vary  with  different  materials  tested,  the  elastic  con¬ 
stants  obtained  may  slightly  differ  from  each  other.  In  the  fol¬ 
lowing  analysis,  the  elastic  constants  of  a  nickel-base  superalloy 
NisAl  single  crystal  given  by  Yang  ( 1985)  are  used.  The  elastic 
compliance  moduli  in  /GPa  with  respect  to  the  crystal  axes 
[001],  [010]  and  [100]  in  the  same  notation  as  Yang’s  are 

Sn  =  Sun  =  1.01 
S12  =  5„22  =  -0.393 
S44  =  45i212  ”  0.848 

The  crystal  orientations  with  respect  to  the  loading  axes 
(fixed  coordinates)  are  represented  by  the  Euler  angles  (0,  </>, 
i//)  as  shown  in  Fig.  3.  The  global  elastic  compliance  in  tensor 
notation  is 


(/>  sin  if/,  -cos  9  cos  4>  smi//  -  sin  <j>  cos  if/,  sin  0  cos  <l> 

</>  sin  t/f,  —cos  0  sin  <f>  sin  ij/  +  cos  </>  cos  if/,  sin  0  sin  <f> 

sin  9  sin  i//  cos  9 


Sijkm  —  ha^jP^k\^my^apky  (23) 

where 

The  inverse  of  Sijkm  is  the  elastic  stiffness  modulus.  Equation 
(23)  was  used  to  obtain  the  elastic  stiffness  matrices  of  single 
crystals  referring  to  the  specimen  axes. 

A  range  of  test  data  of  Ae  Young’s  modulus  and  shear  modu¬ 
lus  of  NisAl  polycrystal  and  its  alloys  are  found  in  the  article 
by  Stoloff  (1989).  The  isotropic  elastic  constants  of  Poisson’s 
ratio,  u  =  0,3,  and  shear  modulus,  ^  -  65  GPa,  are  here  used 
for  our  analysis  of  NisAl  polycryst^s. 

5.2  Physical  Model  and  Numerical  Results.  A  physical 
model  to  study  the  elastic  anisotropy  and  mhomogeneity  of  crystals 
on  the  fatigue  crack  initiation  in  a  polycrystal  is  shown  in  Fig.  4. 
An  aggregate  of  crystals  with  various  crystallographic  orientations 
is  embedded  near  the  free  surface  in  a  semi-infinite  elastic  medium. 
The  cross-sections  of  these  crystals  are  either  whole  or  part  of  a 
hexagon.  Since  the  crystals  at  some  distance  away  firom  the  surface 
crystal,  in  which  a  fatigue  band  is  analyzed,  are  randomly  oriented, 
these  surrounding  crystals  are  assumed  to  be  isotropic  and  homoge¬ 
neous.  The  fatigue  band  is  composed  of  three  thin  slices  P,  Q, 
and  R.  The  orientation  of  surface  crystal  #1  is  determined  by  the 
orientation  of  the  primary  slip  system,  which  is  45  deg  inclining 
to  X2-axis.  The  elastic  stiffness  matrix  of  this  crystal  with  respect 
to  the  loading  axes  is 
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244.33 

118.19 

83.91 

118.19 

244.33 
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83.91 

83.91 

278.61 

-13.71 

41.14 

-27.43 

-41.14 

13.71 

27,43 

.  -6.86 

-6.86 

13.71 

-27.42 

-82.29 

-13.71 

82.29 

27.43 

-13.71 

-54.86 

54.86 

27.43 

153.56 

27.43 

27.43 

27.43 

153,56 

-27.43 

27.43 

-27.43 

222.13. 

The  surface  crystal  with  the  fatigue  band  and  the  surrounding 
crystals  are  shown  in  Fig.  4.  The  fatigue  band  is  divided  into 
parallelogram  elements  and  the  remaining  parts  of  the  crystals 
are  divided  into  trapezoidal  elements.  The  paralellogram  ele¬ 
ments  are  assumed  to  have  linear  inelastic  strains  and  the  trape¬ 
zoidal  elements  to  have  constant  inelastic  strains.  Fig.  5  shows 
the  layout  of  these  . elements  in  the  aggregate. 

The  local  stresses  in  the  surface  crystal  containing  the  fatigue 
band  are  affected  not  only  by  the  applied  load  but  also  by  the 
orientations  of  its  surrounding  grains.  Plastic  strains  can  develop 
both  inside  and  outside  the  fatigue  band  and  the  residual  stresses 
are  affected  by  the  orientations  of  the  surrounding  crystals.  The 
global  elastic  stiffness  matrix  of  ciystals  of  the  orientation  {9, 
<t>,  il/)  is  denoted  as  C*(^,  0, 

The  surrounding  crystals  can  have  different  sets  of  orienta¬ 
tions,  In  the  following,  one  set  is  calculated  to  illustrate  the 
effect  of  anisotropy  and  inhomogeneity  on  the  fatigue  crack 
initiation.  In  this  set,  the  orientation  of  crystals  #2  and  #3  is 
chosen  to  be  (60,  30,  30  deg),  crystals  #4  and  #5  to  be  (0,  0, 
0  deg)  and  crystal  #6  to  be  (30,  30,  60  deg).  The  global  elastic 
stiffness  matrices  of  crystals  with  these  orientations  are  listed 
as  follows 


the  sliding  starts  inside  the  band,  where  the  Schmid  factor  takes 
the  maximum  value.  The  case,  in  which  cr^2  =  391  MPa,  = 
±0.5 MPa  and  =  200MPa,  was  analyzed.  The  surface  plastic 
resolved  shear  strain  vs.  the  number  of  cycles  is  shown  in  Fig. 
7.  If  crystal  anisotropy  is  neglected,  the  fatigue  band  will  not 
develop  in  the  surface  crystal  under  the  applied  load.  Therefore, 
the  threshold  of  fatigue  crack  initiation  is  lowered  with  the 
consideration  of  crystal  anisotropy.  Fig.  8  shows  the  plastic 
resolved  shear  strain  distributions  in  P  at  100  cycles  and  2000 
cycles,  respectively. 

The  cumulative  surface  plastic  shear  strains  versus  the 
number  of  loading  cycles  for  two  cases,  one  with  the  consid¬ 
eration  of  crystal  anisotropy  and  one  without,  were  calcu¬ 
lated  and  plotted  in  Fig.  9.  It  is  seen  that  the  neglect  of  the 
crystal  anisotropy  in  a  polycrystal  may  cause  significant  er¬ 
ror  in  calculating  the  slip  distribution  and  in  estimating  the 
fatigue  crack  initiation. 
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C*(60,  30,  30  deg)  = 
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With  the  method  described,  the  modified  Schmid  factor  Sap 
of  the  primary  slip  system  varying  along  slice  P  is  shown  in 
Fig.  6.  The  Schmid  factor  varies  from  0.42  at  the  free  surface 
to  the  maximum  value  of  0.59.  It  is  seen  that  the  effect  of  this 
anisotropy  on  the  response  to  the  applied  load  makes  the  Schmid 
factor  greater  than  0.5  in  a  large  part  of  the  slice.  The  drop  of 
the  inodified  Schmid  factor  close  to  the  free  surface  seems  to 
be  due  to  the  less  constraint  on  the  deformation  near  the  surface 
than  in  the  interior.  Without  considering  the  crystal  anisotropy, 
the  Schmid  factor  is  uniformly  0.5.  The  increase  of  the  Schmid 
factor  decreases  the  applied  load  required  to  activate  the  plastic 
slip  in  the  band  in  comparison  with  the  isotropic  and  homoge¬ 
neous  case.  The  region  in  which  the  total  resolved  shear  stress 
reaches  the  critical  shear  stress  will  slide  first.  Therefore, 
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